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ABSTRACT 



A spectral analysis of the nearby isolated neutron star RX J1856.5-3754 is presented. Applying the kinetic approach, the distribution 
functions of emitting electrons are derived and the entire spectra is fitted. It is found that waves excited by the cyclotron mechanism, 
come in the radio domain. We confirm that the cyclotron instability is quite efficient, since the estimations show that the time of 
wave-particle interaction is long enough for particles to acquire perpendicular momentum and generate observed radiation. The lack 
of rotational modulation is discussed and the pulsar spin period is estimated. 
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1. Introduction 

The isolated neutron star RX J1856. 5-3754 (hereafter 
RXJ 1856) was discovered by ROSAT as an X-ray source 
dWalter et al. 1996T l. According to the observational evidence 
the emission of RXJ 1856 does not show any significant periodic 
variations and the featureless X-ray spectrum can be fit by 
the Planckian spectrum with a temperature A;r^ ^ 63 + 3eV 
(Burwitz et al. 2003). It has been proposed that the emission of 
this source has a thermal nature and every effort is bent to make 
a model which would well describe the overall spectra. 

First spectral modeling of RXJ 1856 has been presented by 
[Pavlov et al. (1996)] It was shown that the light element (hy- 
drogen or helium) nonmagnetic NS atmosphere models can be 
firmly ruled out, because they over-predict the optical flux by 
a large factor. On the other hand no acceptable fit can be ob- 
tained with iron and standard solar-mixture atmosphere models, 
because the features predicted by these models are not detected 
with high significance. Doppler smearing of the spectral lines 
due to fast rotation of a NS does not wash away completely 
the strongest spectral features ( |Pavlov, Zavlin & Sanwal 20021 
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The similar problems appear 
tized NS atmosphere models 
IZavlin & Pavlov 2002l i. So one has 
classic NS atmosphere models are 
observed X-ray emission of RXJ 1856. 

Different explanation have been considered in papers: 
[Burwitz et al. (2001)HTurolla, Zane & Drake (2004)| It has been 
proposecfl that the star has no atmosphere but a condensed mat- 
ter surface. The mentioned surface might result in a virtually 
featureless Planckian spectrum in the soft X-ray band. But yet 
another problem arose from the fact that the parameters derived 
from X-rays do not fit the optical spectrum with an intensity 6 
times larger, than that of a X-ray emission. This situation led 
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[Pons et al. (2002)[ to introducing the overall spectra by two com- 
ponents. In this model the soft-component of A:r^ ^ 20eV rep- 
resents radiation from relatively cool surface and fit the optical 
data, when the hard component of A:r^ ^ 55eV emitted from 
~ 20% of the NS surface is being responsible for the X-ray emis- 
sion. 

According to [Burwitz et al. (2003) 1 the simplest way to pro- 
duce a time constant flux would be to assume a uniform tem- 
perature distribution across the stellar surface. In this model the 
optical data match the Planckian spectrum for a blackbody tem- 
perature kT^i^ ^ 63eV and a radius ^ 12.3km respectively, 
but the X-ray emissivity is below that of a blackbody by a sub- 
stantial factor. To explain this fact and also fit the overall spectra, 
|Burwitz et al. (2003)| supposed that the radiating surface should 
have a high reflectivity in the X-ray domain. 

So to make these models work, it has been postulated that 
the star has a condensed matter surface. The existence of such 
a medium requkes low temperature and strong magnetic field 
(Lai& Salpeter 1997 ILai 200 il l. The former seems to be ful- 
lilled for both models (Pons et al. 2002, Burwitz et al. 2003), 
when the magnetic field of the star is still unknown. There are 
no convincing theories of existence of such a medium, in which 
it would be possible to confirm the required non-uniform dis- 
tribution of the surface temperature (two-component blackbody 
model). Also it should be shown that such surface has a required 
reflectivity in the X-ray domain, which is one of the main con- 
ditions in the model of [Burwitz et al. (2003)[ Existing models, 
based on an assumption that the emission is of thermal nature, 
have a number of problems. 

We are not about to reject the existing models, but in present 
paper we propose our explanation of emission of RXJ 1856, 
based on well developed models of pulsars. The lack of de- 
tected spin-modulation is easily explained by assumption of a 
nearly aligned rotator (see Fig. [Tjl, in this case the line of sight 
of an observer always lies in the emission cone and the ob- 
server receives continuous radiation. We will show below that 
the presented model does not face with a problem of the ab- 
sence of spectral features in the observed radiation. The lack 
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Fig. 1. O is pulsar spin-axis, ji -axis of magnetic-momentum 
06 « 1°) 

of the spectral lines imposes very stringent constraints on ther- 
mal radiation models. However, a recent paper by Ho et al. 
(2007) explains the observed featureless spectra of RXJ1856, 
based on an assumption that the star has a thin magnetic, par- 
tially ionized atmosphere on the top of the condensed sur- 
face. Though creation of thin hydrogen atmospheres appears to 
be one of important problems for this model. To explain the 
measured spectra we refer to the pulsar emission model pro- 
posed by |Lominadze, Machabeh & Mikhailovskii (1979)| and by 
Machabeli & Usov (1979), where it has been shown that the 
waves excited by the cyclotron instability interact with parti- 
cles, leading to the appearance of pitch-angles, which obviously 
causes synchrotron radiation. 

In this paper, we describe the emission model (§2), derive X- 
ray (§3) and optical (§4) spectra, based on the emission model 
presented in §2, confirm the efficiency of the cyclotron mecha- 
nism (§5) and make conclusions (§6). 

2. Emission model 

It is well known that the distribution function of relativistic 
particles is one dimensional at the pulsar surface, because any 
transverse momenta (p^_) of relativistic electrons are lost during 
a very short time(< 10"^°s) via synchrotron emission in very 
strong B ~ 10"G magnetic fields. The distribution function is 
shown in Fig. |2] (lArons 1981I I. For typical pulsars the plasma 
consists of the following components: the bulk of plasma with an 
average Lorentz-factor y,, ^ 10; a tail on the distribution func- 
tion with Yt - 10"^ and the primary beam with y/, ^10^. Though, 
plasma with an anisotropic distribution function becomes unsta- 
ble that can lead to a wave excitation in the pulsar magneto- 
sphere. The main mechanism of wave generation in plasmas of 
the pulsar magnetosphere is the cyclotron instability. The pro- 
cess can be conveniently described in cylindrical coordinates 
(see Fig. |3]l. The cyclotron resonance condition can be written 
in the form ( [Kazbegi, Machabeh & Mehkidze 1991) : 

Ll)b 

a> - k^V^ - kjcUj, + — = 0, (1) 

where ^ c(l - § - 5^), k^ + kl^ k^, kl = k] + fe^, y, is the 

Lorentz-factor for the resonant particles, m,i is the drift velocity 
of the particles due to curvature of the field lines: 




Fig. 2. Distribution function of a one-dimensional plasma in the 
pulsar magnetosphere. Left corresponds to secondary particles, 
right to the primary beam. 

where p is the radius of curvature of the field lines and oib = 
eB/mc is the cyclotron frequency. The resonant condition for 
transverse (t) waves with the spectrum: 



LJi — kc(l — 6), 
where 

takes the form: 



1 {k±/k^-ujc) \]A 
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This equation describes the wave excitation process by the 
anomalous Doppler eff'ect. During the quasi-linear stage of the 
instability a diffusion of particles arises not only along but also 
across the magnetic field lines. Therefore, plasma particles ac- 
quire transverse momenta (see Eq.l2) and as a result the syn- 
chrotron mechanism is switched on. 

The kinetic equation for the distribution function 
of the resonant particles can be written in the form 
(IMachabeli & Usov 1979IIMachabeh et al. 20021) : 

+ + Gil)/] + ^ ^ [-A (Gx + FJf] = 



+ ^\D\u^ + Duj-]f {pM (6) 
dp\\\ di// dp\\l 

where G is the force responsible for conservation of the adiabatic 
invariant p\lB{r) =const, F- is the radiation deceleration force 
produced by synchrotron emission and Q\\ is the reaction force 
of the curvature radiation. They can be written in the form: 



mc 

G± = Jrl//, 

P 



Gil ■■ 
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where as - 2e w^/Sc and Uc 
will take the form: 



2e2/3p . For the/j, the Eq. (6) 
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the stationary state (df/dt = 0). Then Eq. (12) gives the mean 
value of the angle i^: 



o.04rr 



1/2- 



(15) 



As a result of appearance of the pitch-angles the synchrotron 
emission is generated. 
Using Eq. (15) we get; 
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Fig. 3. Coordinate system (p is a local radius of curvature) 

where = e^/mc^, tf/Q is the mean value of the pitch angle and 
\Eic[^ is the density of electric energy in the waves. 

The resultant spectrum of synchrotron emission is defined by 
the shape of distribution function of the resonant particles, which 
in turn is the solution of Eq. (10). To define the total flux emitted 
by the resonant particles we use the power of a single electron 
multiphed by the distribution function of emitting electrons and 
integrate the derived expression over the energy distribution (see 
[Ginzburg 1981] ). 



3. X-ray spectrum 

We suppose that the measured X-ray spectrum is the result of 
synchrotron emission of primary beam electrons. Let us assume 
that before the cyclotron instability arises, the energy distribu- 
tion in the beam has a shape: 



rib 



exp 



(r - jh) 



(11) 



where yb - 10^, yj- - 10 - is the half width of the dis- 
tribution function and Ub - BjPce is the density of pri- 
mary beam electrons, equal to the Goldreich-Julian density 
dOoldreich & Juhan 1969l l. 

The wave excitation leads to a redistribution process of the 
particles via the quasi-linear diffiision. Let us compare relative 
values of all the three terms in the right-hand side of Eq. (10). 
Here we consider that y/^i/'o » 1 (the validity of this assumption 
will be estimated after obtaining the expression for (/tq )• In this 
case the equation for the diffusion across the magnetic field has 
the form (Malov & Machabeh 2002): 



df 



bo 



1 



dt mcyji// dtff 



1 d L^^Jfb. 
m^c^yjiff dt// \ dt// 



where 



D 



\Ek\\ 



2 1 mc^ribybC 
2 cj 



(12) 
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(14) 



The transversal quasi-linear diffusion increases the pitch- 
angle, whereas forces Ij and ~P resist this process, leading to 



where P is the pulsar rotation period. We suppose that the period 
of RXJ1856 is f ~ Is (this is a typical period for the pulsar with 
the age~ lO^yr). Using the parameter values indicated above, we 
will take that the quantity (16) is of order of 10. 

Now let us compare the second and the third terms in the 
right-hand side of Eq. (10): 



27T^reiJ/a\Ek\^/y,. 
acyj 



10' 



11/4 

rb 



g2p3/4y5^,9/2- 



(17) 



We consider that the magnetic field is dipolar Let us assume 
that for the primary-beam electrons the resonance condition is 
fulfilled at the distance r ^ lO'cm then for the magnetic field 
at the same distance, we will have B ~ lO^G. For the parame- 
ter values pointed above, it turns out that, the third term of the 
Eq.(lO) is by six orders of magnitude less than the second one. 
Then one gets: 



(18) 



For ytt// « 10'", a magnetic field inhomogeneity does not 
affect the process of wave excitation. The equation which de- 
scribes the cyclotron noise level, in this case has the form: 



d\Ejl_ 
dt 



= 2r\E,\%, 



where F - is the growth rate of instability. 



F = : 
and 
F = , 
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where the resonant frequency 
( [Kazbegi, Machabeli & Mehkidze 199 l| l 



is defined 



Sy,' 



(19) 



(20) 



(21) 



(22) 



In this case we can use Eq. (21) for the growth rate and after 
substituting Eq. (19) in (18), we will take: 



d_ 

dt 



f\\ 2 2^yb^\Ek\ 

^ m cot, u'- oyb 



= 0. 



(23) 



From this equation it is easy to find the distribution function: 



f\\b -f\\i,o + ^Tfc- 



(24) 
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We suppose that the right slope of the measured X-ray spectrum 
is the result of synchrotron emission of the primary-beam elec- 
trons, with the distribution function (24), see Fig.|4] 

In |Malov & Machabeh (2002)| it has been shown that it 
is possible to form the stationary distribution. If we assume 
dfu/dt = 0, then from Eq. (18) one finds: 



-4 



(25) 



We consider that after the quasi-linear evolution stage of the in- 
stability, by achieving the stationary state, the radiation density 
is significant and the self-absorption effects begin to play the 
main role. Synchrotron self-absorption, redistributes the emis- 
sion spectrum and in the domain of relatively low frequencies, 
we have (Zh elezniakov 1977 J : 

I(v) oc v^/^. (26) 

This emission corresponds to energy interval (0.15-0.26 keV) of 
the measured X-ray spectrum (see Fig.|4|. 

The frequency of the maximum of measured X-ray spectrum 
corresponds to the intersection point of theoretical curves, fitted 
with the observed data (v ^ 6 ■ 10"'Hz). On the other hand, the 
frequency of the maximum of synchrotron spectrum of a single 
electron is ( Pacholczyk et al. 1970[ l: 



eB 2 
Vm - 0.07 — sin\py,. 
mc 



(27) 



Using Eq. (15) for \p «i 1, we find that the frequency of max- 
imum of the spectrum turns out to be not dependent on magnetic 
field. Consequently, we can estimate from Eq. (27) the rotation 
period of RXJ1856. Estimations show that P ~ Is, i.e. the prim- 
itive conjecture proves to be true. The same can be found about 
two other assumptions done at the beginning of our computa- 
tions, since i/^oi ~ lO"-' <s; 1 and jb^ot ~ lO-'. 

4. Optical spectrum 

Now let us consider the observed optical emission of RXJ1856. 
We suppose that this emission is related with electrons of a tail, 
with an average Lorentz-factor y, ^ 10"*. Let us assume that the 
initial distribution of particles of the tail has the following form 
(see Fig.|2]i: 



(28) 



As in the previous case, one needs to estimate the contribu- 
tion of different terms in the right-hand side of Eq. (10). We sup- 
pose that the condition y,iAo » 1 is fulfilled and we can use Eq. 
(15) for the mean value of angle t/f in this case too. For values of 
parameters indicated above we take that, the first term exceeds 
the second one by five orders of magnitude. Now let us compare 
the first and the third terms: 



2n^reipQ\Ek?-ly, 

.2,,2 



1.3- 10^- 



3,3/2" 



(29) 



We suppose that for the electrons of the tail, the cyclotron res- 
onance condition is fulfilled at the distance r ^ 3 ■ lO^cm and 
for the value of the magnetic field we take B ~ \G. Finally, after 
substituting values in Eq. (29) we take that the third term is less 
than the first one by an order of magnitude. 

Consequently, instead of Eq. (10), we will have: 

In^^or, d (\E, 



dt 



dy\ y 



0. 



(30) 



In this case we can use Eq. (20) for the growth rate of instability 
and from Eqs. (19) and (30) one gets: 



d_ 
Jt 



' n^e^cDByT , d \Ek\ 



2\ 



dy y^ 



0. 



(31) 



From this equation we can find the expression for \Eit\^ at the 
final moment of the quasi-linear stage. If we take into account 
that the density of the particles behaves as l/P (r is the distance 
from the pulsar), one can neglect /j| in comparison with/j| g and 
from Eq. (31) one finds: 



\Ek 



2 m^r'doj^p^ 



iTT^e^ajgyTi^Qy' 



(32) 



After achieving the stationary mode, the distribution func- 
tion of particles of the tail can be found from Eq. (30) if one sets 
dfu/dt - 0, then we will take: 



.^1 



\Ek\-h. 



Using Eq. (32) for the \Ek\ we find: 

.,2 



(33) 



(34) 



f\\, r 

The resultant theoretical spectrum is well matched with the 
measured one. For parameters used above, we find that indeed 

r,iAo, ~ 10^ » 1 (-Ao, - 3 ■ 10-'). 

The frequency of the maximum in synchrotron emission of a 
single electron is defined by Eq. (27). For the parameter values 
used above, we will take v„j(Optic) ^ 2 ■ lO'^Hz. This quantity 
comes in the domain that corresponds the observed optical data. 

5. The effectiveness of the cyclotron mechanism 

Now let us estimate the frequencies of original waves, excited 
during the cyclotron resonance. Using the same parameters, one 
can see from Eq. (22) that the frequencies corresponding to the 
beam and tail components are v a; 3 ■ lO'^Hz and v » 5 ■ 10*^ Hz 
respectively, i.e. during the cyclotron resonance in both cases the 
radio waves are excited. Though this waves are generated earlier 
than X-ray and optical emission and therefore the radio emission 
might pass by line of sight. In this case the radio emission will 
not get to observer That is one of possible explanations why the 
radio emission is not detected from this object. 

For effective generation of waves it is essential that time, 
during which the particles give energy to waves should be 
more than l/F. Generated radio waves propagate practically 
in straight lines, whereas the field line of the dipole mag- 
netic field deviates from its initial direction and the an- 
gle 6 - k\\lky. grows {6 is the angle between the wave 
line and the line of dipole magnetic field). On the other 
hand, the resonance condition (5) imposes limitations on 9 ^ 
max{ V^, — ) (Kazbegi, Machabeli & Melikidze 1991|l, i.e. parti- 



cles can resonate with the waves propagating in a limited range 
of angles. Obviously the Eq. (5) will be fulfilled until then 
( [Lyutikov, Blandford (fe Machabeli 1999| l: 



c 



(35) 



where cF ' is the growth length and 9p is the length of the wave- 
particle interaction. For the beam particles from Eq. (35) it fol- 
lows that p > 3 ■ 10''cm. As the cyclotron instability arises at 
distances r ^ lO^cm for the beam electrons, this result means 
that time of wave interaction with the resonant particles is quite 
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Fig. 4. The Chandra LETG X-ray spectramofRXJl 856.5-3754 
fitted with the model. 

enough for particles to acquire the pith-angles (/'o (Eq. 1 5), which 
automatically leads to the generation of synchrotron emission. 

The total energy available for the conversion into pulsar 
emission is of the order of the energy of primary beam particles 
flowing along the open filed lines of the pulsar magnetosphere: 

E - nh,nRl^ybmc\ (36) 

where nb^ is the Goldreich- Julian density at the pulsar surface 
and Rpc is the radius of the polar cap. The estimations show 
£■ ~ 5 . lO^^erg/s, which is quite enough to explain the observed 
luminosity of RXJ1856. 

6. Conclusions 

1 . In the present paper we assume that the emission of RXJ 1 856 
is generated by the synchrotron mechanism. The main rea- 
son of wave generation in outer parts of the pulsar mag- 
netosphere is the cyclotron instability. During the quasi- 
linear stage of the instability a diffusion of particles arises 
as along, also across the magnetic field. Plasma particles ac- 
quire pitch angles and begin to rotate along the Larmor or- 
bits. Synchrotron emission is the result of the appearance of 
pitch angles. 

2. The measured X-ray spectrum is the result of synchrotron 
emission of primary-beam electrons, generated at distances 
r ^ lO^cm, when the optical emission is related with sec- 
ondary plasma particles, particularly the tail electrons. The 
optical spectrum is produced by synchrotron emission of the 
tail electrons, generated at distances r ^ 3 • lO^cm. The pre- 
dictable characteristic frequencies Vm{X - ray) ^ 6 ■ lO'^Hz 
and v,n{Optic) ^ 2 - 10"*Hz come in the same domains as the 
measured spectra. 

3. The observed optical emission reveals an intensity six times 
larger than the X-ray one. The X-ray emission is generated 
earlier in contrast to optical emission. Thus, an observer may 
receive only part of the radiation emitted in the X-ray do- 
main. We suggest that this must be the reason of lower inten- 
sity of the X-ray spectrum. 

4. In view of the fact that the aligned rotator model is used, the 
radio emission covers a large distance in the pulsar magne- 
tosphere. So there is a high probabihty for it to come in the 



cyclotron damping range w - k^Vy, - kxUx - ^ = 0. In this 
case the radio emission wiU not reach an observer. Another 
explanation of the lack of radio emission from this object is 
that the X-ray and optical emission, resulted from the quasi- 
linear diffusion, occur in a region further away from that of 
the radio emission and it miss our line of sight. Nevertheless 
the detection of radio emission from RXJ1856 would be 
strong argument to our model. 
5. We give estimations of effectiveness of the cyclotron mecha- 
nism, which in this case means the fulfillment of the follow- 
ing condition p > 3 • lO'cm. As the instability develops at 
distances r ^ lO^cm, then follows that the excited waves lie 
in the resonant region sufliciently long time, which is quite 
enough for particles to acquire perpendicular momentum and 
generate radiation. 
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